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ABSTRACT 

We present results from a study of the globular cluster luminosity function (GCLF) in a sample of 89 early- 
type galaxies observed as part of the ACS Virgo Cluster Survey. Using a Gaussian parametrization of the GCLF, 
we find a highly significant correlation between the GCLF dispersion, a, and the galaxy luminosity, Mga-iX, in 
the sense that the GC systems in fainter galaxies have narrower luminosity functions. The GCLF dispersions 
in the Milky Way and M3 1 are fully consistent with this trend, implying that the correlation between sigma 
and galaxy luminosity is more fundamental than older suggestions that GCLF shape is a function of galaxy 
Hubble type. We show that the cr-Mfi gai relation results from a bonafide narrowing of the distribution of 
(logarithmic) cluster masses in fainter galaxies. We further show that this behavior is mirrored by a steepening 
of the GC mass function for relatively high masses, A^>3xlO^A40,a mass regime in which the shape of the 
GCLF is not strongly affected by dynamical evolution over a Hubble time. We argue that this trend arises from 
variations in initial conditions and requires explanation by theories of cluster formation. Finally, we confirm 
that in bright galaxies, the GCLF "turns over" at the canonical mass scale of A^to — 2 x IQr' M.q. However, we 
find that A^jo scatters to lower values (« l-2x 10^ A^©) in galaxies fainter than M^ gai > -18.5, an important 
consideration if the GCLF is to be used as a distance indicator for dwarf ellipticals. 

Subject headings: galaxies: elUptical and lenticular, cD — galaxies: star clusters — globular clusters: general 



m 
o 
o 

Oh 

6 



1. INTRODUCTION 

The luminosity function of globular clusters (GCs) repre- 
sents one of the most remarkable features of these stellar sys- 
tems. The distribution of GC magnitudes, commonly referred 
to as the GC luminosity function (GCLF), shows a turnover, 
or peak, at My — -7.5 mag, corresponding to a mass of 
A( ~ 2 X IQt'AAq. Observations have shown that this turnover 
is nearly invariant across and within galaxies, prompting its 
widespread use as a distance indicator (see, e.g., Harris 2001). 
Accounting for this nearly universal mass scale remains an 
open problem for theories of GC formation and evolution. It 
follows that establishing whether or not the GCLF as a whole 
is universal — i.e., whether its overall form depends on host 
galaxy properties — can help guide and constrain theories for 
the formation and evolution of galaxies and GC systems. 

In this Letter, we present results from a study of the GCLFs 
of 89 early-type galaxies observed by HST as part of the ACS 
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Virgo Cluster Survey (ACSVCS; Cote et al. 2004). We find 
the clearest evidence to date for a correlation between the 
width (i.e., Gaussian dispersion) of the GCLF and the lumi- 
nosity of the host galaxy; we also show that there is some 
downward scatter in the mass scale of the GCLF turnover 
in galaxies fainter than Mb.^hx > -18.5. Focusing on the ob- 
served steepening of the GCLF at the bright (high-mass) end 
in the faint galaxies, we argue that this behavior was prob- 
ably imprinted at the time of GC formation. A more de- 
tailed discussion of the whole GCLF, including the faint (low- 
mass) end and the role that long-term dynamical evolution 
plays in that regime, is deferred to a subsequent paper (Jordan 
et al. 2006, hereafter J06). That paper presents our data in full 
and gives details of our analysis techniques, including model- 
ing of the GCLFs with a new, non-Gaussian, physically moti- 
vated fitting function. 

2. OBSERVATIONS AND ANALYSIS 

One hundred early-type galaxies in the Virgo cluster were 
observed in the ACSVCS (Cote et al. 2004). Each galaxy was 
imaged for 750 s in the F475W bandpass (~ Sloan g) and for 
1210 s in F850LP (~ Sloan z). Reductions were performed as 
described in Jordan et al. (2004). In what follows, we use g 
and z as shorthand to refer to the F475W and F850LP filters. 

One of the main scientific objectives of the ACSVCS is the 
study of GC systems, and thus we have developed methods to: 
( 1 ) discard foreground stars and background galaxies from the 
totality of observed sources around each target galaxy in the 
survey; and (2) estimate the level of residual fore- and back- 
ground contamination in the remaining sources designated as 
candidate GCs. These procedures are described and illus- 
trated by Peng et al. (2006a; their §2.2 and Figure 1), and 
discussed in detail in the GCLF context in J06. In the latter 
paper, we also examine the effects of using alternate selec- 
tion criteria to define GC samples, and show that the results 
presented here are fully robust against such subtleties. 

Of the 100 galaxies in the ACSVCS, we restrict our anal- 
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ysis to those that have more than five GCs, as estimated by 
subtracting the total number of expected contaminants from 
the full list of GC candidates for each galaxy. We additionally 
eliminate two galaxies for which we were unable to obtain 
useful measurements of the GCLF parameters. This leaves a 
final sample of 89 galaxies which are studied here, and in J06. 

Also as part of the ACSVCS, we have measured the dis- 
tances to 84 of our target galaxies using the method of sur- 
face brightness fluctuations (SBF; Mei et al. 2006). We use 
these SBF distances to transform the observed GC and galaxy 
magnitudes into absolute ones whenever possible. For those 
galaxies lacking an SBF distance, we adopt the mean dis- 
tance modulus to the Virgo cluster: ((m-M)o) = 31.09 mag, 
or {D) = 16.5 Mpc (see Mei et al. 2005, 2006). 

We use an approach similar to that of Seeker & Harris 
(1993) to characterize the GCLFs: parametric models are 
fitted to the observed luminosity functions via a maximum- 
likelihood method that takes into account photometric errors, 
incompleteness, and the luminosity function of contaminants. 
Full technical details are given in J06, where we consider two 
parametric models for the GCLF. The first, on which this pa- 
per will focus, is the standard Gaussian distribution, 

dN/dz = N,o,{2TTa^y^'^ exp[-(z-/i,)V2a?] . (1) 

The second is a simple analytical modification of a Schechter 
(1976) function designed to account for the effects of clus- 
ter evaporation (two-body relaxation) on a GC mass function 
that is assumed to have initially resembled that of the young 
clusters forming today in local mergers and starbursts. Full 
details on these two models are given in J06, where we fit 
each of them to the separate g- and z-band GCLFs of our 89 
program galaxies. In this Letter we present only the results of 
Gaussian fits to the z-band GCLFs. 

3. RESULTS 

Figure ^ shows our main result: GCLFs are narrower in 
lower-luminosity galaxies. The straight line in this plot of 
Gaussian dispersion against absolute galaxy magnitude shows 
the least-squares fit 

cr, = (1.12±0.01)-(0.093± 0.006)(MB,gai+20) . (2) 

It has been reported before that the GCLFs in lower- 
luminosity galaxies tend to show somewhat lower dispersions 
(e.g., Kundu & Whitmore 2001). However, the size and ho- 
mogeneity of the ACSVCS dataset make this the most con- 
vincing demonstration to date of a continuous trend in GCLF 
shape over a range of > 400 in galaxy luminosity. Monte 
Carlo simulations and alternate constructions of GCLF sam- 
ples show that the observed decrease in dispersion is not an 
artifact of small-number statistics in the faint galaxies (J06). 

Past investigations have pointed to a dependence of the 
GCLF dispersion on the Hubble type of the GC host galax- 
ies (e.g., Harris 1991). Figure [l] includes datapoints at the 
location of the bulge magnitude and GCLF dispersion of the 
Milky Way (large star) and M31 (large triangle). Since both 
systems fall comfortably on the relation defined by our data 
for early-type galaxies, we conclude that the underlying fun- 
damental correlation is one between a and Mg gai, rather than 
between a and Hubble type. 

A natural question at this point is whether the observed 
trend in GCLF dispersion with galaxy magnitude implies a 
similar trend in the GC mass function. This is not a foregone 
conclusion, for the following reason. GC systems are known 
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Fig. 1. — Gaussian dispersion, a-, versus galaxy, Mg gal. for the z-band 
GCLFs of 89 ACSVCS galaxies. Tlie GCLF width varies systematically, 
being narrower in fainter galaxies. The two anomalously high points at 
Ms.gA =-21.2 and -19.9 correspond to the galaxies VCC 798 and VCC 2095, 
both of which have large excesses of faint, diffuse clusters (Peng et al. 2006b). 
The large star is plotted at the spheroid luminosity (de Vaucouleurs & Pence 
1978) and GCLF dispersion (Hairis 2001) of the Milky Way. The large tii- 
angle marks the bulge luminosity (Kent 1989) and GCLF dispersion (Hanis 
2001)ofM31. 



to have systematically redder and broader (or more strongly 
bimodal) color distributions in brighter galaxies than in faint 
ones (see, e.g., Peng et al. 2006a). Equivalently, GCs in gi- 
ant galaxies are more metal-rich on average, and have larger 
dispersions in [Fe/H], than those in low-mass dwarfs. Since 
cluster mass-to-light ratios, T, are functions of [Fe/H] in gen- 
eral, it is conceivable that the average GC T could change 
systematically in going from bright galaxies to fainter ones, 
and that the spread of T values within a single GC system 
could also vary systematically as a function of galaxy mag- 
nitude. The possibility then exists that narrower GCLFs for 
faint galaxies might result from these systematics in T com- 
bined with a more nearly invariant spread in GC masses. We 
can show easily, however, that this is not the case. 

The systematics in T vs. [Fe/H] just mentioned are also a 
function of wavelength. In bluer filters, such as B, V, or g, 
mass-to-light ratios of old stellar systems do change signifi- 
cantly (increasing by factors of two or more) in going from 
cluster metallicities [Fe/H] < -2 to [Fe/H] = 0, typical of 
GCs. But at the much redder wavelengths of our z-band data 
(Apivot — 9055 A; Sirianni et al. 2005), this strong metallic- 
ity dependence almost completely disappears. We have used 
the PEGASE population-synthesis model of Fioc & Rocca- 
Volmerange (1997) to compute T- as a function of metallicity 
for clusters with a Kennicutt (1983) stellar IMF and various 
fixed ages r. For r = 13 Gyr, we find that ~ 1.6 AAqLq 
at an extreme [Fe/H] = -2.3, decreasing to a minimum of 
Tj ~ 1.5 AAqLq at [Fe/H] ~ -0.7, and then increasing 
sUghtly to T, = 1 .7 X at [Fe /H] = 0. In other words, we 
always have T.. « 1 .6 ± 0. 1 for any of the globular clusters in 
any of our sample galaxies — no matter how red or blue the 
clusters are, or how broad or narrow the GC color/metallicity 
distribution. Comparably small ranges of T, result if younger 
GC ages or different reasonable stellar IMFs are assumed. 
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Fig. 2. — (Top) GCLF turnover magnitude (absolute mean fi-) versus 
galaxy magnitude, A/fiuai, from Gaussian fits to 89 z-band GCLFs in the 
ACSVCS. The outlying points at Mg.g^i = -21.2 and -19.9 are VCC 798 
and VCC 2095, which have large excesses of faint, diffuse clusters (Peng et 
al. 2006b). The star and triangle show fi- values for the Milky Way and M31, 
estimated from their V-band peaks (Harris 2001) by applying an average 
(V — z) color estimated from the PEGASE population-synthesis code (Fioc 
& Rocca-Volmerange 1997). (Bottom) Turnover mass Mjq corresponding 
to the fitted fj,-, obtained by applying an average GC T; computed for each 
galaxy using the PEGASE model. Typical errorbars on A4to as a function 
of galaxy magnitude are indicated. 

The effect of variations in mass-to-light ratio on the width 
of the GCLF at NIR wavelengths is therefore completely neg- 
ligible. From Figure^ we have that (T(log L.) = a./2.5 > 0.2 
in our galaxies, whereas the discussion above implies that 
Tj;(max)/T,(min) ~ 1.13 for our GCs, such that the dis- 
persion of mass-to-light ratios in any one system is always 
cr(log Tj.) < 0.055 at an absolute maximum. The intrinsic dis- 
persion of logarithmic GC masses, (T(log A4) = [a^(log L,)- 
cr^(log T,)]'/^, is thus never more than different from 

the observed cr(log L.). We conclude, unavoidably, that the 
steady decrease of by more than 50% from the brightest 
giants to the faintest dwarfs in Figure ^is an accurate reflec- 
tion of just such a trend in the intrinsic GC mass distributions. 

We now turn our attention to the GCLF turnover magnitude. 
The upper panel of Figure |2 shows the mean GC absolute 
magnitude ^~ from the Gaussian fits to our GCLFs, versus 
host galaxy Mg.gai- The horizontal line in this plot is drawn 
at a level typical of galaxies brighter than Mfi gai < -18.5: 
/i; = -8.4. Given a typical ~ 1.5 AAqLq in these galax- 
ies (for GC ages 13 Gyr and an average [Fe/H] « -1), this 
corresponds to a cluster mass scale of AAjo — 2.2 x 10^ .Mq- 
Estimates of the z-band GCLF turnovers in the Milky Way 
and M31 are shown by the large star and large triangle, as 
in Figure ^ In the lower panel of Figure |2l we plot the 
turnover masses A^to obtained from the fitted /ij. using PE- 
GASE model mass-to-light ratios. As we have discussed, z- 
band luminosities are very good proxies for total GC masses, 
so this graph is essentially a mirror image of the one above it. 

Figure |2| shows that there is no strong or systematic varia- 
tion in /i, or A^TO to match that seen for cr, (Figure^. Nev- 



ertheless, there is a clear tendency for the GCLF turnovers 
of galaxies fainter than Mg.gai > -18.5 to scatter to somewhat 
fainter (less massive) values than is typical of the bright gi- 
ants. The difference in mass is a factor of « 1.5 on average, 
but it ranges apparently randomly, from a factor of 1 (i.e., no 
difference) up to factors slightly greater than 2 in some cases. 
Note that there is a healthy mix of E and SO or dE and dSO 
galaxies at all magnitudes in our ACSVCS sample (see Table 
1 of Cote et al. 2004). We find no tendency for any particular 
Hubble type to scatter preferentially away from /i- = -8.4 or 
Mm = 2.2 X IO-'^Mq in Figure|2l 

The lower mean value for A^to at faint Ms.gai clearly 
can impact the use of the GCLF as a standard candle for 
dwarf galaxies. On the other hand, the effect is wavelength- 
dependent. Publicly available codes such as PEGASE can be 
used easily to show that in bluer bandpasses such as g (or the 
closely related V, which is more standard for such studies), 
the slight decrease we find for the average GC turnover mass 
in fainter galaxies is balanced by a comparable decrease in 
the typical GC mass-to-light ratio (because of the lower clus- 
ter metallicities), so that the mean turnover magnitude does 
not vary as in the z band. We have also confirmed this directly 
from our own ACSVCS data. In J06, we obtain plots anal- 
ogous to Figures [2 and 12] from fits to the g-band GCLFs of 
our galaxies. The results fully support all of our conclusions 
here. It is particularly worth noting that we find a relation 
identical to equation (0 for the dependence of g-band GCLF 
dispersion on parent galaxy luminosity. 

4. DISCUSSION 

An obvious question prompted by Figure [0 is whether the 
correlation between a~ and M^.gai was established at the time 
of cluster formation or built up afterwards as GCLFs were 
modified by the dynamical destruction of GCs over a Hubble 
time. We favor the first interpretation. 

Star clusters can be destroyed over Gyr timescales as a re- 
sult of mass loss driven by stellar evolution, dynamical fric- 
tion, gravitational shocks, and internal two-body relaxation 
(evaporation) — processes that have been studied in detail by 
several groups. Recent discussions, centered specifically on 
how these affect the GCLF, can be found in Fall & Zhang 
(2001) and Vesperini (2000, 2001). FaU & Zhang in particu- 
lar show that, while stellar evolution and gravitational shocks 
certainly deplete the total number of GCs in a galaxy, they do 
not significantly alter the overall shape of the GCLF. Evapora- 
tion, on the other hand, can change the shape of the GCLF, but 
significantly so only for cluster masses Ai < 2-3 x 10^ Mq, 
i.e., below the typical GCLF turnover mass. 

In the theoretical treatments of Fall & Zhang, Vesperini, and 
many others, the evaporation rate is independent of cluster 
mass, which ultimately drives the low-mass side of the GCLF 
to a universal shape: a simple exponential dN/dz oc lO"""^' 
for the number of GCs per unit magnitude fainter than the 
turnover (equivalent to a flat distribution for the number of 
GCs per unit linear luminosity or mass). But fitting a Gaus- 
sian model to the GCLF, as we have done here, tacitly as- 
sumes that the distribution is symmetric. The results in Figure 
Hmight therefore seem to imply that both the bright side and 
the faint side of the GCLF become progressively steeper in 
fainter galaxies. However, various observational uncertainties 
make it difficult to determine precisely the form of the faintest 
tail of the GCLF. Thus, in J06 we show that good fits to our 
GCLFs can also be obtained using an alternate model with 
a universal exponential shape at magnitudes fainter than the 
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Fig. 3. — Slope of the power law that best fits our z-band GCLFdata, fi-, for 
masses 3 X lO' <^ (M/Mq) ^ 2 x 10*, plotted against host galaxy absolute 
magnitude, Mbm^x ■ The stai' and triangle show /3 values for the Milky Way 
and M3 1 respectively, measured in the same mass regime using the data from 
Harris (1996) and Reed et al. (1994) assuming a V-band mass-to-light ratio 
M/Lv = 2. The bright side of the GCLF is steeper in fainter galaxies. 

turnover — and that the downward scatter in A^to for faint 
galaxies persists in such a model and so is not an artifact of 
any assumed Gaussian symmetry. Here we concern ourselves 
only with the brighter half of the GCLF, which is observation- 
ally better defined. 

We have performed maximum-likelihood fits of exponen- 
tial models dN jdz oc 10°-^*''""'^~ (corresponding to power-law 
mass distributions, dN jdM, cx MT^-) to the GCLFs at abso- 
lute magnitudes -8.7 > z > -10.8 (cluster masses ~ 3 x 10^- 
2 X 10^ A^o) in 66 of our galaxies. Such distributions accu- 
rately describe the bright sides of giant-galaxy GCLFs (Harris 
& Pudritz 1994; Larsen et al. 2001), and with A ~ 2 they also 
give good matches to the mass functions of young star clusters 
in nearby mergers and starbursts (Zhang & Fall 1999). 

Figure|3lshows the results of this exercise. There is a clear 
steepening in the power-law exponent, from /S^ ~ 1 .8 in bright 



galaxies to /3, ~ 3 in the faintest systems. However the faint 
side of the GCLF behaves in detail, the bright side alone sug- 
gests that smaller galaxies were unable to form very massive 
clusters in the same relative proportions as giant galaxies. 

A potential complication here is dynamical friction. A clus- 
ter of mass hA on an orbit of radius r in a galaxy with circu- 
lar speed Vc will spiral in to the center of the galaxy on a 
timescale Tdf oc A^"'r^Vc (Binney & Tremaine 1987). In the 
Milky Way and larger galaxies, Tdf > 13 Gyr for all but the 
very most massive clusters at small radii, and thus dynam- 
ical friction does not significantly affect their GCLFs (e.g.. 
Fall & Zhang 2001). In dwarfs with low Vc, however, Tdf can 
be interestingly short for smaller GCs at larger r — suggest- 
ing, perhaps, that the process might significantly deplete the 
bright side of the GCLF in small galaxies and contribute to 
the type of trend seen in Figure|3] However, Vesperini (2000, 
2001) has modeled the GCLF evolution over a Hubble time 
in galaxies with a wide range of mass, and his results strongly 
suggest that dynamical friction does not suffice to explain our 
observations. In particular, the widths of the Gaussian GCLFs 
in his models do not decrease, even in dwarf galaxies, to any- 
where near the extent seen in the data. Thus, any significant 
galaxy-to-galaxy variations in the shape of the GCLF above 
the turnover mass probably reflect initial conditions (see J06 
for further discussion). 

In summary, the gradual narrowing of the GCLF as a func- 
tion of galaxy luminosity — or the steepening of the mass dis- 
tribution above the classic turnover point — presents a new 
constraint for theories of GC formation and evolution. In 
our view, it is the cluster formation process in particular that 
is likely to be most relevant to the observed behavior at the 
high-mass end of the GCLF. Exactly what factors might lead 
to more massive galaxies forming massive clusters in greater 
relative numbers, is an open question of some interest. 
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